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We investigated the mutual induction effects between the d-wave and the s-wave components 
of order parameters due to superconducting fluctuation above the critical temperatures and 
calculated its contributions to paraconductivity and excess Hall conductivity based on the 
two-component stochastic TDGL equation, ft is shown that the coupling of two components 
increases paraconductivity while it decreases excess Hall conductivity compared to the cases 
when each component fluctuates independently. We also found the singular behavior in the 
paraconductivity and the excess Hall conductivity dependence on the coupling parameter which 
is consistent with the natural restriction among the coefficients of gradient terms. 
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§1. Introduction 

Recently many experimental results indicate the mixed pairing symmetry of Cooper pairs in high- 
Tc cuprates which consists of d- and s-wave- like components.BS Theoretically, it is shown that 
certain classes of microscopic models lead to dominant d-wave pairing with subdominant s-wave 
pairing.^ Macroscopic properties of superconductors of these mixed symmetry can be studied by 
means of the Ginzburg-Landau (GL) thoery for two-component order parameter. Phenomenologi- 
cally, JoyntS first introduced the GL free energy of two-component superconductors based on the 
group theoretical consideration.^ Microscopic derivations of the corresponding two-component GL 
equations have been performed by several authors.!' 

Due to those macroscopic arguments, it is shown that, in tetragonal lattice, the coexistence of the 
two components is not expected in the bulk, but the subdominant order parameter is induced by 
spatial variations because two components are coupled through "mixed gradient" term; this effect 
has been intensively investigated under the circumstances such as vortices, Si'S'S^ impurities,!-* 
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currentsEi) etc., thereby their anisotropic features being clarified. Such inductions have been also 
studied from more microscopic standpoints. El' Nonetheless there is so far no considerations 
of the effect beyond these static cases. There should be dynamical situations where the order 
parameter varies spatially so that the subdominant component is induced. We now point out this 
possibility, which is driven by the superconducting fluctuaion above critical temperature, and give 
an estimation of the effect to the paraconductivity and the excess Hall conductivity. It should be 
stressed that, within our considerations, the origin of the induction is superconducting fluctuation 
itself, unlike the recent studies treating vortex states based on two-component TDGL equation, 
where the induction is attributed to intrinsic structures of vortices. 

The conduction process which contributes to the enhancement of electrical conductivity due to 
superconducting fluctuaion is first pointed out by Aslamasov and Larkin(AL). Besides, there is 
a prediction about the behavior of Hall conductivity resulting from AL processJlB' though it is not 
yet clear whether it accounts for the Hall anomaly .0 AL process can also be described within the 
framework of Time-Dependent Ginzburg-Landau(TDGL) equation and several formulations have 
been don^ll^ill' based on the stochastic TDGL equation. 

In this work, based on the stochastic two-component TDGL equation, we show our calculation for 
the paraconductivity and the excess Hall conductivity and discuss the importance of the induction 
effect of the subdominant component through the superconducting fiuctuation. Since we limit 
the considerations to the case of weak magnetic field, we will make the simple linear or bilinear 
expansion of the response current in terms of the external electric and magnetic field. 

The characteristic features of the induction effect we will show is firstly that no anisotropy 
arises in spite of the anisotropy in GL free energy and secondly that the induction reduces the 
excess Hall conductivity while it enhances the paraconductivity. The results will show the singular 
behavior of these transport coefficients in increment of the coupling of two components through 
the spatial variation. This is direct refiection of the stability condition of the normal state against 
the superconduting fluctuation. 

In Sec.|, we first give the formulations of fluctuation conductivity tensor based on the two- 
component TDGL equation as a stochastic differential equation. In Sec.^, we show the numerical 
results and see how the coupling of two components affects them. In Sec.^, we summarize the 
results and give some discussions. In this paper we set h = kB = c = 1. 

§2. Formulation for the stochastic TDGL equation 

In this section we show the formulation for the paraconductivity and excess Hall conductivity 
based on the stochastic two-component TDGL equation. In the following, electric and magnetic 
field is expressed as scalar and vector potential (j){r),A{r), respectively. The charge of a Cooper 
pair is e* = — 2e < 0. We introduce a differential operator 77 = — iV — e*A{r). 
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We start with the two-component TDGL equation 



7.(^ + ie0jd = -(^ + a.U '-^^s + f, (2.1a) 



e e 
jf(r) = — (d*nd + C.C.) + — (s*ns + C.C.) 

rrid rus 

e* e* 
-X — {(fn^s + s*n^4 + c.c.) + y — {(fllys + s*nyd + c.c.) (2.2) 



which is based on Joynt's two-component GL free energylf and also can be derived from such lattice 
models as t-J model.El' d and s mean superconducting order parameters of d-wave and s-wave 
component, respectively. Both GL coefficients and Og are positive because the system is above 
the transition temperatures. We added the white noise source terms /i,/2 as representations of 
thermodynamical fluctuations, which make the mean square expectation value of order parameter 
finite even above the transition temperature. It is the main characteristics of the TDGL equation 
that the two components are coupled via so-called "mixed gradient" terms. The main purpose 
of this paper is to clarify the effects of this term on the electric transport properties, and this 
effects are represented as the coefficient m~^. Because of this anisotropic term we should have set 
electrical field as E = {E^, Ey, 0) for more generality, but it will be found later that this anisotropy 
is cancelled with that of current and no anisotropy arises in the fluctuation conductivity tensor. 
For this reason we set (j){r) = —Ex, A{r) = Bxy. 

First, we rewrite the TDGL equation( |2.1| ) into the simpler form. Keeping up to bilinear terms 
in external fields, Fourier transformed TDGL equation becomes 

\-iuJ'4){k,u:) = - {r}k + lVk^\ V'(fc,w) + f{k,uj), (2.3) 



where 



Vk = e*{E-\B^-^kyWy) (2.4) 



(2.5) 
(2.6) 



% = '^Wx + -^Wy + a, (2.7) 
W.= [ \ \], (2.8) 
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Wy=[ :j, (2.9) 

-(::)■ 

Then we can formally express the solution of the TDGL equation: 

il:{k, Lo) = G{k, io)f{k, io) (2.12) 

where 



G{k,Lo)=(^iju; + r]k + 7Vk-^^ . (2.13) 



d_ 

Here we assume the correlation function between white noises as 

{fik',u')fHk,u;)) = aS (2.14) 

where the coefficient a is 2 x 2 matrix and should be determined from the linearized two-component 
GL free energy 

Fk = ilj{k)^mHk). (2.15) 

ip{k) means the order parameter of wave vector k at the specific time. Corresponding distribution 
function is ^ 

Zk= f e-^'=/^ddidd2dsids2 = TT-^' (2-16) 
J det T]k 

where indices 1 and 2 mean real and imaginary part, respectively. This distribution function gives 
rise to 

{mm% = -T { f^;^ 'l^^ 1 InZ, = Tv^\ (2.17) 



d 


1 d 






1 d 


d 


2 drjv 


dria 



{■ ■ •)o means thermal average in the absence of external fields. The stochastic TDGL equation ( |2.lD 
also realizes this equilibrium values, provided we set 

a = r(7 + 7t). (2.18) 

In fact, from eg. ( 2. 14 ) with eq. (|2.18| ), we obtain 

(V(fe)V'(fc)t)o = ^(V'(fc,u;)V^(fc,a;))o 

= Y,Go{k,u){f{k,Lo)p{k,uj))Gl{k,u;) 



TVk' (2-19) 
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where Go(fc,w) = {i^io + %)^"^ stands for the Green's function of TDGL equation( |2.3| ) in the 
absence of external fields. 

Next we perform expansions with respect to external fields V^. Resolvent G{k,ij) is expanded as 
follows: 

G{k,u;) = Go{k,uj) + Gi{k,uj) + G2{k,u;) + ■ ■ ■ . (2.20) 

Note that, in the above expansion, each term acts on any function h{k) as a differential operator 
except Go{k,uj): 

Gi{k,uj)h{k) = -Go(fc,u;)7Vfc — [Go(fc,a;)Mfc)], (2.21) 



G2{k,u:)h{k) = GQ{k,u:)^Vk^ 



Goik,ojhVk^[Go{k,oj)h{k)] 



(2.22) 



In contrast, G^(fc,Ci;) acts on functions at the left side. 

On the other hand, the expectation value of supercurrent(2.2) is written as the thermal averaged 
form as 

iJir)) = e* E E E E ^'mik) - e*i)]e^(^'-^)- (2.23) 

k k' Lu' 

where 

(m, k 771.7 fc\ ~ 

' \ ,fc = (-fe,,A;,), (2.24) 
k rrig k j 

A=[ \ \ ,A = (-^.,^,), (2.25) 

and 

R{k,io;k',Lo') = {tlj{k' ,uj')tlj\k,Lo)). (2.26) 

Corresponding to the resolvent expansion ( |2 . 2C| ) , the correlation function R{k,uj; k' ,uj') is expanded 
as follows: 

R{k,uj;k',uj') = G{k' ,J){f{k' ,J)f\k,uj))G\k,uj) 

= Ro{k, io; k', to') + Ri {k, u; k', cj') + i?2(fc, w; k' , to') + ■ ■ ■ (2.27) 

where 

Ro{k,oj;k\J) = {2TT)^Go{k' ,uj')a5{k - k')S{uj - uj')Gl{k,Lo), (2.28) 



Ri{k,uj;k',uj') = (27r)^[Gi(fc',cj>5(fc - k')6{uj - u')Gl{k,uj) 
+Go{k',uj')a6{k - k')5{oj - J)G\{k,uj)], 



(2.29) 



R2{k,iv;k',u;') = {2TT)^[G2{k' ,uj')a6{k - k')6{uj - uj')Gl{k,u;) 

+Gi{k',uj')a6{k - k')5{oj - J)G\{k,uj) 
+Go{k', uj')a5{k - k')5{oj - oj')G\{k, uj)]. 



(2.30) 
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Now we can calculate supercurrent from the correlation functions. The expressions for the su- 
per current in each order are as follows: 



{Joir)) = E E E E t^[^o(fe, u;; k', J)k{k)]S''-^)\ (2.31) 

k k' u)' 

(j 1 (^)) = e* E E E E [^1 {k, u:- k', J)k{k) - Ro{k, k\ J)e* A]S^' -^^^ (2.32) 

k k' ^ uj' 

{h{r)) =e*J2T.T.T. tr[^2(fe, u;; k' , J)k{k) - R, (k, u;; k\ J)e* A]S^' (2.33) 
k k' u}' 

Here we consider a thin film with thickness d, so that we calculate the summation as integral with 
respect to kx and ky, and extract only the contributions from kz = 0. Thus we have the expression 
for the paraconductivity from ( p.32| ) 

CT.. = ^EEt^[i-^^'(4^-G'oW^.Go7Go + h.c.) + gS7GSw^.Go}(7 + 7^)] (2.34) 

u> k 

and for the excess Hall conductivity from ( p. 33] ) 
ie*^TB 

<Jyx = E E tr[[{A;'M^j,(l - 2klGoWx)GoWxGo{WyGoj + jGoWy) - h.c.} 

^ k 

-{klklWyGoWxGoiWyGoWxGo^ + jGoW^GoWy) - h.c.} 

+{klkl{WyGlWxGlWyGoWxGo-/ - h.c.}]Go(7 + 7^)4] (2-35) 

where Go means Go{k,uj). 

From the above expression we find that both axx and ayx are even functions of m^^ and therefore 
we obtain ayy = axx, CTxy = —<7yx\ no anisotropy arises in the fluctuation conductivity tensor in 
contrast to recent theoretical predictions about vortex state below the critical temperatureJll'llll' 

Since the above expression for conductivity tensor is complicated, it may be helpful to compare 
it with the former studies in the case of = 0. In this case the coupling of two components 
disappears and the problem must reduce to that in the conventional single component super- 
conductor, which has been studied Aslamasov and Larkin(AL)i) and Fukuyama, Ebisawa and 
Tsuzuki(FET).ll3' Then all matrices become diagonal and actually one can confirm 

(Jxx = E E E(^i + Ii)[-^li^i'^li9*i9l + c.c.) + m-^j^gfgi] 

i=d,s <^ k 

= y flplL (2.36) 



and 



ie*^TB 

oyx = E EE(7i + 7;)[2/c2m-37,gf5* 

i=d,s ^ k 

-6klk'^m-^-figfg* + klk^m-^-figfgf - c.c] 
_ eW^^ (2.37, 



6 



where we introduced single component Green's functions 

1 e , , 

gi = : ,Vi = 7, ^ "i- (2-38) 

and 

7ii = Re7i, 7,2 = Im7i (2.39) 
These results correspond to AL and FET, respectively. 



Fig. 1. Temperature dependence of (a)paraconductivity axx and (b)excess Hall conductivity ayx for the coupling 
parameter m^^ = O(solid line),m^^ — 3(long dashed line),m^^ = 4 (short dashed line) where Td — 20, Ts = 
16, m^^ — 3, mj^ = 6, 7d = 12 + 4i, 7^ = 4 — lOi. Here we incorporated the temperature dependence by setting the 
GL parameters Ud = T — Td and as = T — Ts- 



§3. Numerical results 

In this section we show the numerical results of paraconductivity ( p . 34 ) and excess Hall 
conductivity(2.35). In Fig.l we show the plot of them as functions of temperature for various 
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coupling parameter m~ . Here the critical temperature of the d-wave component is chosen to be 
higher than that of s-wave, and as a result, both Gxx and ayx diverges at even in the presence 
of the coupling of the two component order parameters. It S66II1S tllclt (7 XX is enhanced but Oyx is 
surpressed due to the coupling of the two component order parameters. In order to see the coupling 
effects more clearly, we show the above results as functions of in Fig. 2. In the light of GL free 
energy ( p^ ), detW^^ > is required for the stability of the system against the spatial variation 



of order parameters. This condition is reflected in the singularities at = ity rus in the 
fluctuation conductivity tensor. 



Fig. 2. Coupling parameter dependence of (a)paraconductivity a^x and (b)excess Hall conductivity ayx at fixed 
temperature T = 22 with the same set of parameters as in Fig.l. 



§4. Summary and Discussions 

In this work we gave the formulation for the paraconductivity and the excess Hall conductivity 
in the presence of the coupling of two-component order parameter in the limit of weak external 
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fields. The phenomenological formulations based on the stochastic TDGL equation have been done 



considered to be their extension. We also numerically calculated the dependences on temperature 
and coupling parameter with the typical parameters and found that the coupling of two components 
enhances paraconductivity but reduces excess Hall conductivity. Such tendency does not change 
even for the other set of parameters. 
We paid attention to the stability condition about GL coefficients and saw the singular behavior 



of conductivity tensor at critical value ~ -^y ' though it could have not been attained 

if one had resorted to perturbative expansions with respect to m~^. 

It is well known that excess Hall conductivity vanishes identically in the absence of imaginary 
part of order parameter relaxation time within the conventional framework of TDGL theory. 0) On 
the other hand, several recent studies investigating vortex states based on two-component TDGL 
equationJli'0' have pointed out that the Hall effect is also caused by the anisotropy of crystal 
lattice even if the relaxation times are real. However, our calculation in the fluctuation regime 
resulted in vanishing Hall conductivity in the case of Im7rf = Im7s = 0, and moreover, we found 
that no anisotropy arose in the conductivity tensor even for finite value of m~^. 

Finally we mention an aspect of multiple transition in unconventional superconductors from 
the viewpoint of superconducting fluctuation. If the transition temperatures are close to one 
another in the presence of the coupling between components, the induction effect in superconducting 
fluctuation would gain its signiflcance. In this context, the recent experiments through microwave 
surface impedancei^ which suggest multi-component superconductivity can present an interesting 
system if those components correspond to order parameters of different symmetry. To compare with 
such experiments, we need to carry out the extention so as to include the frequency dependence in 
paraconductivity, which will be given elsewhere. 
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